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^Nj Abstract. Violent relaxation - the protocluster dynamical response to the expulsion of its resid- 
ual star forming gas - is a short albeit crucial episode in the evolution of star clusters and star 
O cluster systems. Because it is heavily driven by cluster formation and environmental conditions, 
. it is a potentially highly rewarding phase in terms of probing star formation and galaxy evolu- 
tion. In this contribution I review how cluster formation and environmental conditions affect the 
shape of the young cluster mass function and the relation between the present star formation 
rate of galaxies and the mass of their young most massive cluster. 
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1. Introduction 

Compact gas-embedded clusters are observed to be the formation sites of a significant 
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^ . fraction of stars in the focal Universe. Following star formation, feedback processes - 
[ ranging from protostellar outflows, photoionization-driven gas overpressure, stellar winds 
^ • to, eventually, supernovae - expel the gas not yet converted into stars. That is, the 
' efficiency with which a cluster parent core turns its gas into stars is always smaller than 
^— I . 100 per cent. Not only does gas expulsion terminate star formation, it also leaves the 
^ ' newly-formed cluster out-of-equilibrium, thereby forcing it into violent relaxation. In an 
attempt to regain a state of equilibrium, the cluster expands, loses a fraction of its stars 
_ or, even, undergo complete disruption. These are the cluster infant weight-loss and cluster 
• infant mortality, respectively. 

■ The dynamical response of clusters to gas expulsion implies that the cluster age distri- 
I bution, the cluster mass function, the ratio between the mass in clusters and the mass in 
stars, as well as the distribution function of cluster half-mass radii all constitute promis- 
ing diagnostic tools of star cluster formation conditions. Baumgardt & Kroupa (2007) 
^ ' have built and made available to the community an A'^-body model grid describing the 
• ^ ' temporal evolution of cluster mass-loss and cluster spatial expansion for a wide range of 
local star formation efhciencies (SFE; i.e. the mass fraction of dense molecular gas that 
. cluster- forming cores turn into stars), gas-expulsion time-scales and external tidal field 
■ " " ' impacts. To date, this is the most comprehensive grid of model clusters through violent 
relaxation. Combined to the ever increasing quality of cluster observational data-sets 
(including the much needed cluster half-mass radius, see e.g. Scheepmaker et al. 2007), 
the time is now ripe to decipher observed properties of young star clusters so as to probe 
their formation and environmental conditions. 

The outline is as follows. Following a brief reminder of violent relaxation chief pa- 
rameters in Section [2l I highlight how formation and environmental conditions affect 
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the shape of the cluster mass function in Section [31 I also present some preliminary re- 
sults about ongoing work regarding the relation between the luminosity of the brightest 
young cluster and the present star formation rate (SFR) of spiral galaxies, namely, the 
SFR-M^"^''*''' relation (Section H]). 

2. Gas expulsion: prime parameters 

In what follows, the mass of a cluster is defined as the stellar mass enclosed within its 
tidal radius. We assume that the impact of the external tidal field arises solely from the 
smooth gravitational potential of a Milky- Way-type galaxy, i.e., at this stage, we neglect 
the potential impact on the early evolution of a star cluster of its natal giant molecular 
cloud. We adopt the Jacobi radius of a cluster as an estimate of its tidal radius (see 
Baumgardt et al., in press, their eq. 1). 

In the course of violent relaxation, the mass TOci of a star cluster evolves rapidly and 
must therefore be defined as an instantaneous quantity. It obeys: 

nicl = Fb{t/Tcross, £,TGE^p/Tcross, Th/n) X E X 711^ . (2.1) 

In this equation, rric is the mass of the cluster-progenitor molecular core, e is the 
local SFE and is the stellar mass fraction of the initially gas-embedded cluster still 
contained within the cluster tidal radius at time t. When using the term '(gas-)embedded 
cluster mass', this refers to the stellar mass med = sx nic. In all the work presented here 
the local SFE is assumed to be uniform among all cores of a given simulation (i.e. no 
core-mass dependent SFE). The bound fraction Fi, depends on the time t/rcross elapsed 
since gas expulsion and on the gas expulsion time-scale tq Exp /''cross (both expressed in 
unit of the molecular core crossing-time), on the local SFE e and on the impact of the 
external tidal field. The latter is quantified as the ratio between the half-mass radius r/j 
and the tidal radius of the gas-embedded cluster (see below for details). The older the 
cluster and/or the lower the SFE and/or the quicker gas expulsion and/or the stronger 
the host galaxy tidal field, the lower the bound fraction Fi,. Building on the motion of 
a supershell propagating through the cluster-forming core and collecting the residual 
star-forming gas, Parmentier et al. (2008) derived an expression for the gas-expulsion 
time-scale tgexp/tctoss (their eq. 6; see also fig. 3 in Parmentier & Fritze 2009). Because 
a cluster expands following gas expulsion, the mass fraction of stars it loses depends on 
how compact the gas-embedded cluster is compared to its tidal radius rt. The impact of 
the external tidal field is therefore accounted for by the ratio Vh/rt. As will be shown 
in Section [31 the cluster-forming core mass-radius relation is therefore of paramount 
importance to the outcome of violent relaxation since it influences both the depth of the 
core potential well, hence the gas-expulsion time-scale, and the sensitivity to the external 
tidal field of the exposed cluster (see also Parmentier 2009). 

3. Young cluster mass functions 

In this section, I show how specific combinations of the core mass-radius relation, local 
SFE and external tidal field lead to cluster mass functions markedly different from the 
core mass function and, therefore, to mass-dependent cluster infant mortality. Cluster- 
forming cores are spatially resolved in the nearby regions of our Galaxy only and the 

f We explicitly assume that stars have had sufficient time to come into virial equilibrium with 
the gas potential (i.e. gas expulsion occurs after a few crossing-times) and hence that the local 
SFE matches closely the effective SFE (see Goodwin (2009) for a discussion). 
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Figure 1. Time-evolution of the cluster mass function when cluster-forming core radii are 
constant (rc = 1 pc) and the local SFE is 20 per cent. The symbol- free dotted (black) line is the 
power-law embedded-cluster mass function ('ECIMF'). The lines with crosses, open circles and 
filled squares (blue, pink and red, respectively) show the cluster mass functions integrated over 
the age ranges quoted in the key. 



slope of the core mass-radius relation remains therefore ill-determined (see Parmentier, 
Kroupa & Baumgardt, subm. for a discussion). I consider two extreme cases in turn: 
constant core radii and constant core surface densities. 

Constant cluster-forming core radii and mass-varying gas expulsion time-scale. 

Under the assumption of constant radii (rc = 1 pc), a sequence of increasing core mass 
equates with a sequence of increasing densities and deepening potential wells. As a result, 
high-mass cores undergo slower gas expulsion and retain a larger fraction of bound stars 
than their low-mass counterparts. Under the assumption of a weak tidal field, Parmentier 
et al. (2008) showed that the shape of the post-violent-relaxation cluster mass function is 
strongly sensitive to the local SFE (see their figs. 2 and 3). While a local SFE e ~ 0.4 is 
conducive to a cluster mass function mirroring the core mass function (mass-independent 
infant mortality and weight-loss), a local SFE e c± 0.2 results in the survival of the - 
initially - highest-mass embedded-clusters only. For an exposed cluster to survive so 
low a local SFE, gas must be expelled adiabatically, a condition which is met by high- 
mass cores and their deep potential well. Only those produce bound clusters at the end 
of violent relaxation, rendering cluster infant mortality heavily mass-dependent. As a 
result, a power-law core mass function gives rise to a bell-shaped cluster mass function. 

An aspect which Parmentier et al. (2008) did not investigate is the time-evolution 
of the power-law mass function of embedded clusters in the bell-shaped cluster mass 
function. The topic of evolving secularly (i.e. when violent relaxation effects have faded 
away) a power-law cluster mass function in a bell-shape is extensively covered in the 
literature. Starting from a power-law at an age of ~50-100 Myr, internal 2-body relaxation 
and external tidal stripping preferentially removes low-mass clusters, thereby carving a 
turnover in the mass function. As time goes by, this turnover moves towards higher 
cluster mass and smaller mass function amplitude. However, the evolution of relevance 
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Figure 2. Time-evolution of the bound fraction of stars for an SFE of 20 per cent and 
two distinct regimes of cluster gas expulsion: explosive (i.e. rcExp << Tcross) and adiabatic 
(tgbxp ~ 3 Tcross). Under the assumption of constant core radii, high-mass clusters have a higher 
density, which results in slower gas expulsion, shorter crossing-time and - once exposed - faster 
evolution than their low-mass counterparts. 

here is violent relaxation instead of secular evolution and, as we shall see, the cluster 
mass function evolves differently. 

This is the perfect example to illustrate that the evolutionary stage of a cluster is not 
determined by its age in an absolute sense. As for the gas expulsion time-scale, what 
matters is the age expressed in unit of a core crossing-time (see eq. 12. ip . This implies 
that, keeping all other parameters alike, clusters born in cores of higher densities evolve - 
once exposed - at a rate faster than those formed out of low-density cores. Figure [1] shows 
the embedded-cluster mass function (i.e. the core mass function shifted by —0.7 in logm 
to account for the constant 20 per cent local SFE) and snapshots of the cluster mass 
function at ages of ~ 15, 30 and 45Myr. The early cluster mass function (age ~ 15Myr, 
dotted line with cross-signs) shows prominent substructures, with the cluster mass still 
little affected by infant weight-loss in the low- and high-mass regimes {rrici < 3, 000 Mq 
and rrici > 2 x 10® M0, respectively). Their origin arises from how the bound fraction Fj, 
evolves with time and as a function of gas-expulsion time-scale. In Fig. [2l the solid and 
dotted lines show the Fb time-evolutions for explosive (tqexp << Across) and adiabatic 
(here tqexp = 3 Tcross) gas expulsions, respectively. A massive embedded cluster (say, 
rrieci = 2 X 10® M©) expels its residual gas on an adiabatic time-scale and hence retains 
a bound core of stars in spite of the low SFE of 20 per cent. Also, by virtue of the high 
density of its natal core, a young age of 15Myr equates with several hundreds of core 
crossing-times. Accordingly, cluster violent relaxation is over by this young age. This 
is why most of the bell-shape, which arises from clusters initially more massive than 
TOecJ ^ 10® Mq, is already carved by an age of 15Myr. 

In contrast, low-mass clusters (e.g. nieci = 200 M©) have low core densities, long core 
crossing-times and the same age of 15 Myr equates with a few Tdoss only, a stage at which 
the instantaneous bound fraction Ff, is still close to unity. Consequently, the low-mass 
regime of the cluster mass function at an age of 15 Myr sticks to the embedded-cluster 
mass function. As clusters age to 50 Myr, the cluster mass function in the high-mass 
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Figure 3. Dependence of the ratio between tiie lialf-mass radius and tlie tidal radius rt of 
the embedded cluster on its stellar mass rrif-ci for constant surface density cores. Owing to their 
large spatial extent, high mass clusters show a higher sensitivity to the external tidal field and 
present a greater likelihood of disruption following gas expulsion. 

regime does not evolve anymore and the low-mass regime of the mass function disap- 
pears as low-mass objects experience infant-weight loss until their eventual dissolution. 

Constant cluster-forming core surface densities and mass-varying tidal field impact. 

Following Parmentier et al. (2008), molecular cores of constant surface density have an 
(almost) constant gas-expulsion time-scale. Yet, the infant mortality and infant weight- 
loss rates of their embedded clusters are markedly mass-dependent owing to the tidal 
field impact. We insist that this effect takes place even if all cores are located within 
the same limited region of a galaxy within which the external tidal field does not vary 
markedly. The ratio between the initial half-mass radius and initial tidal radius of these 
embedded clusters scales as rh/rt oc rn'J'^ I'rfiJ'^ (x rnj^ and, therefore, high-mass clus- 
ters have a greater sensitivity to the external tidal field since the volume of space in 
which they can freely expand is smaller than for low-mass clusters when compared to 
the initial cluster size. This is illustrated in Fig. [3] where the core mass-radius relation is 

~ 0.008 (mc/1 M©)^/^ PC0, the local SFE is e = 0.4 and the galactocentric distance is 
^gai = 5kpc in a Milky- Way-like potential. As a result, a power-law core mass function 
evolves in a post-violent-relaxation cluster mass function depleted in high-mass objects 
(see Fig. HI). 

What does it all mean ? 

The key-point to retain from this section is that the shape of the cluster mass function, 
the mass-radius relation of cluster-forming cores, their local SFE and the external tidal 
field are intimately entwined issues. That most observed mass functions of young clusters 
are power-laws with a spectral index —2 (see Parmentier & Gilmore (2007) and references 
therein; but see also Anders et al. (2007) for a possible turnover in the cluster mass 

f Core mass-radius relations adopted in this contribution are taken from Parmentier & Fritze 
(2009). An update on these relations will be presented in Parmentier, Kroupa & Baumgardt 
(subm.) 
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Figure 4. Constant surface density cores: when more massive objects are more vulnerable. 
The dotted (black) line is the core mass function and the solid (red) line is the cluster mass 
function at an age of ~ 50Myr. As can be understood from Fig. |3l cluster infant- mortality acts 
dominantly in the high-mass regime. Model parameters are the same as in Fig. |3] 

function of the Antennae merger NGC 4038/39) therefore tells us something about cluster 
formation conditions which are common to different types of cluster environments. The 
high sensitivity to the external tidal field of high-mass cores of constant surface density 
may be a hint that, contrary to theoretical expectations, the index /? of the mass-radius 
relation Vc oc of cluster-forming cores is shallower than (3 = 1/2. More simulations 
are required to settle this point definitively. In particular, greater attention must be paid 
to cluster radius-related quantities, e.g. the distribution function of young cluster half- 
mass radii and the reason why clusters lack a clear mass vs. half-mass radius correlation 
(Larsen 2004). We will investigate these aspects in a forthcoming paper. 



4. Brightest young star clusters in galaxies 

Star clusters are chronometers of major star formation events in their host galaxies 
(Kroupa 2002; KotuUa, Fritze & Anders 2008). When the SFR of a galaxy increases, not 
only does the galaxy form a larger number of gas-embedded star clusters, the cluster 
mass function is also sampled up to a higher cluster mass, partly as a result of the 
size-of-sample effect. Time spans in the history of galaxies characterized by vigorous 
star formation therefore show up prominently in cluster age histograms for two reasons: 
(i) the number of clusters is larger and (ii) these more numerous clusters also include 
more massive ones which experience longer dissolution time-scales and which are therefore 
more likely to be observed. 

Observationally, this effect is best-known as a correlation between the present SFR of 
galaxies and the absolute visual magnitude ^^'^ahtest i^j^gjj. brightest cluster (Larsen 
2002; Weidner, Kroupa & Larsen 2004; Bastian 2008). Building on synthetic cluster 
populations, Bastian (2008) shows that i/ galaxies form all their stars in gas-embedded 
clusters, then a universal mass fraction of 8 per cent of clusters survive the transition 
from their gas-embedded stage to being exposed. 

That the ratio between the mass in clusters and the mass in stars depends sensitively 
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Figure 5. Observed relation between tiie absolute visual magnitude of the brightest cluster of 
spiral galaxies and their present SFR (filled triangles). The dotted (pink) line is the relation 
between the SFR and the mass of the most massive core for the adopted Schechter core mass 
function (a = —2, mhr = 1.2i?6MQ). Solid (red) and dashed (blue) lines depict an iso-SFE model 
and its 1 — CT limits. We caution that, at the exploratory stage of this model, the evolutionary 
time span of the model (lOMyr) is narrower than that of the observations (> 100 Myr), thereby 
preventing us from drawing conclusions about the best local SFE estimate (see text for details) 



on the mean local SFE (Parmentier & Fritze 2009) prompted Parmentier, Kroupa & 
Baumgardt (subm.) to investigate what cluster formation conditions are the main drivers 
of the SFR - My"^^*''"* relation. The absolute magnitude of a cluster depends on its 
age-dependent mass-to-light ratio, on the stellar mass of its parent gas-embedded cluster 
and on the amount of infant weight-loss this has experienced. Section [3] demonstrates 
that cluster mass-related quantities are heavily shaped by cluster formation conditions. 
We have therefore divided our study of how the SFR — M^^s^test j.gjg^^jQjj responds 
to input parameter variations in two main parts. Firstly, we consider the limited age 
range of 10 Myr over which we explicitly assume the constancy of the integrated cluster 
mass-to-light ratio, that is, we investigate how the mass of the most massive cluster 
responds to cluster formation conditions, specifically the core mass-function, core mass- 
radius relation, local SFE and strength of an external tidal field (Parmentier, Kroupa & 
Baumgardt, subm.). In a second forthcoming paper, we will extend the cluster age range 
to 100 Myr and assess how mass-to-light ratio variations, hence the age of the observed 
brightest cluster, contribute to shape the SFR ~ M^^shtest ^^i^^^q^^ fQj. ^j^g gj.g^ ^^^.^ 
of the study, our findings are as follow. The vertical location of the SFR — ]\/[^^9htest 
relation depends on the external tidal field and local SFE. Stronger tidal fields and/or 
smaller SFE are conducive to smaller mass/brightness for the most massive/brightest 
clusters. As for the shape of the SFR — ]\/f^'^9htest j-gjg^^jQj-j^ jg dictated mostly by the 
core mass function (see Fig. [5]) and the core mass-radius relation. For constant surface 
density cores, the SFR — ]\/[^^9htest ,-gjg^^jQj-^ gg^g shallower when the SFR increases due 
to the greater sensitivity to the external tidal field of high-mass embedded clusters and 
hence the larger amount of infant-loss they experience. In contrast, for constant core 
radii, it gets steeper towards higher SFR as, in this case, higher mass cores experience 
slower gas expulsion and thus retain a larger bound fraction of stars. 
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In the case illustrated in Fig. [Sj cluster-forming cores have a constant volume density 
{rc = 0.026 (mc/1 Mq)^/'^ pc), the local SFE is 35 per cent and the synthetic cluster 
population is located at a galactocentric distance of 3 kpc in a Milky- Way-like potential. 
The core mass function is a Schechter function with a spectral index —2 and its break- 
mass TTifjr has been adjusted so that, following star formation and infant weight-loss, 
the break-mass of the cluster mass function is about ~ 2 x 10^ Mq as observed for 
spiral galaxies (Larsen 2009). The solid and dotted lines are the model and its 1 — a 
limits, respectively. The filled triangles are the data-points from Larsen (2002) for spiral 
galaxies. Figure [5] suggests that an iso-SFE model provides a good fit to the ensemble 
of the data, provided that other conditions ~ core mass-radius relation (its slope and 
its normalization), external tidal field - do not vary significantly from one spiral galaxy 
to another. We emphasize once again that the simulation time span (lOMyr) is shorter 
than the age range of clusters shown in this plot (> 100 Myr). Since the model vertical 
location is also cluster-age dependent through both the amount of infant weight-loss and 
the cluster mass-to-light ratio, the limited time span of our simulations prevents us from 
drawing firm conclusions about the best local SFE estimate. Having now a firm handle 
of how cluster formation and environmental conditions affect the relation between the 
present SFR and the mass of the young most massive cluster, we are ready to perform 
simulations on a grander scale encompassing ~ 100 Myr of dynamical evolution and in 
which cluster mass-to- light ratio variations will be taken into account. This will give us a 
lower limit on the local SFE in gas-embedded clusters of spiral galaxies, the actual value 
depending on the importance of the distributed mode of star formation. 
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